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We compare two Monte Carlo implementations of resummation schemes for the description 
of parton evolution at small values of Bjorken x. One of them is based on the Balitsky-Fadin- 
Kuraev-Lipatov (BFKL) evolution equation and generates fully differential parton distributions in 
momentum space making use of reggeized gluons. The other one is based on the Catani-Ciafaloni- 
Fiorani-Marchesini (CCFM) partonic kernel where QCD coherence effects are introduced. It has 
been argued that both approaches agree with each other in the x — > limit. We show that this is 
qq ' not the case for azimuthal angle dependent quantities since at high energies the BFKL approach is 

dominated by its zero conformal spin component while the CCFM gluon Green function receives 
' contributions from all conformal spins even at very small x, 

o ' 

1 Introduction 

An important challenge in the phenomenology of Quantum Chromodynamics (QCD) is to under- 
stand what are the dominant effective degrees of freedom underlying the strong interaction at high 
energies. In the limit where the center-of-mass energy y/s in a collision is much larger than any of 
the relevant mass scales the Balitsky-Fadin-Kuraev-Lipatov (BFKL) approach |T}0 appears to be 
a very useful framework to describe the scattering. In its original formulation this approach is based 
on the exchange of "reggeized" gluons in the i-channel. The interaction among them takes place via 
a gauge invariant (reggeized gluon)-(reggeized gluon)-(gluon) vertex. This simple effective struc- 
ture stems from the dominance of the so-called multi-Regge kinematics where gluon cascades are 
only ordered in longitudinal components. 

The simplicity of the final integral equation should not shadow the strong self-consistency of 
the full BFKL program where tight bootstrap conditions linking the reggeization of the gluon with 
the pomeron wavefunction, dominant in diffractive interactions, are fulfilled. The regime of ap- 
plicability of the leading order BFKL approach to phenomenology is limited since it should be 
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valid in a certain window of center-of-mass energies and perturbative scales. To extend its range 
of applicability one should include either higher order corrections, going beyond the 'multi-Regge" 
kinematics, include non-linear corrections responsible for the restauration of unitarity at very small 
aQ or, as we are going to discuss in this paper, to include a more global treatment of collinear regions 
in phase space using the Catani-Ciafaloni-Fiorani-Marchesini (CCFM) equation, which provides a 
good matching from BFKL to the x — > 1 regime, at least, as we are going to show in this work, as 
long as anomalous dimensions and k t -diffusion properties are concerned. 

In this letter we first give a brief introduction to the BFKL and CCFM approaches, in section [2] 
We then describe our Monte Carlo implementations in sections[3]and[U where we also compare both 
resummation schemes in terms of the detailed exclusive information we can obtain from running the 
Monte Carlo codes. In particular, we show that the UV/IR symmetry in the k t -diffusion present in 
the BFKL evolution is broken at lower energies in the CCFM equation, to then be slowly restored 
as x tends to 0. Our main result appears when we study the azimuthal angle dependence of the 
gluon Green function in both approaches. We show that the CCFM equation generates a stronger 
azimuthal angle dependence than the BFKL equation. To stress this point we investigate the Fourier 
components of the solution in this angular sector. We show that the projection on the zero Fourier 
component, or zero conformal spin in the BFKL context, has the highest growth with energy in 
both schemes. The differences appear in the non-zero Fourier components, or non-zero conformal 
spins, which in the BFKL description fall with energy but in the CCFM solution increase. We also 
discuss the different collinear behaviour of the solutions at different energies. Finally we present our 
conclusions and scope for further investigations. 



2 Multi-Regge kinematics versus QCD coherence in small x fi- 
nal states 

Coherence effects are already present in Quantum Electrodynamics where they are responsible for 
the suppression of soft bremsstrahlung from electron-positron pairs. In QCD processes such as 
g —> qq where a soft gluon is emitted from one of the fermionic lines with an angle larger than the 
angle of emission in the original qq pair probes the total colour charge of this pair. This charge is the 
same as that of the parent gluon and the final radiation takes place as if the soft gluon was emitted 
from it. This is what we call a colour coherence effect and leads to the angular ordering of sequential 
gluon emissions. In Deep Inelastic Scattering (DIS) we can focus on the (i — l)th emitted gluon 
with energy Ei—i from the proton. A secondary gluon radiated from it with a fraction (1 — zi) of its 
energy and a transverse momentum q± has the opening angle 

(1) 



(1 - Zi)Ei-x 1 

with Zi = Ei/Ei-i and \9%\ <C 1. Colour coherence leads to angular ordering with increasing 
opening angles towards the hard scale, in this case the virtuality of the photon Q 2 , and we have 

9 i+ i > Oi, or 

flH-l >M (2) 



J- ~ Zi+l 

which in the limit Zi, Zj+i <C 1 is equivalent to > Ziqi. 



In this work we consider equivalent variables the center-of-mass energy s, the DIS Bjorken x and the rapidity Y linking 
them by Y ~ In s ~ In l/x. 
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In Refs. the BFKL equation for the unintegrated gluon density was written in a form 

suitable for the study of exclusive quantities: 

Uk) = /°(fc) + a s / -| / -z"A R (z, k)Q(q - n)f u (q + k), (3) 
J f| ! Jo z 

where u> is the Mellin-conjugate variable to Bjorken x, \l plays the role of a collinear cutoff, q is the 
transverse momentum of the emitted gluon, and the reggeized gluon propagator is 

/ 1 fc 2 \ 

A fl (z i ,fc l ) = cxp -a s ln — ln-| , (4) 
V z i M / 

with fcj = \ki\, and as = asN c /ir. Solving the equation by iteration real gluon emissions are 
generated with the corresponding virtual corrections summed to all orders. Since f u is inclusive and 
IR finite it corresponds to a sum over all final states with the /i-dependence cancelling between the 
real and virtual contributions in the final result. 

The DIS structure function is calculated integrating over all /i 2 < qf < Q 2 and reads 

oo „Q 2 r ,2 

^ BFKL (Q,^e(Q-/i) + E / II-% d *— (5) 



where i corresponds to the number of gluon emissions, which in the leading order approximation 
coincides with the number of iterations of the kernel. In Ref. J3] it was pointed out that coherence 
effects should significantly modify each contribution to the final sum with a fixed number of emit- 
ted gluons, r, whilst preserving the sum. Therefore, care must be taken to account properly for 
coherence in the calculation of associated distributions. 

Modifying the BFKL formalism to account for coherence JUQ, i 7 J FKL becomes 

fS cfm (Q^) = e(Q-M) 

oo „Q 2 r <2 

+ IT — T dz * — ZiA{zi,qi,ki)@(qi —Zi-iqi-i), (6) 

where A(zi, q$, ki) is not a reggeized gluon propagator anymore but stills plays the role of a no- 
emission factor and takes the CCFM form 



A(zi,qi,ki) = exp 



1 , k 



2 

-as In — In 

Zi Ziqf 



h> qi- (7) 



For the first emission q^z^, = /u is chosen. This collinear cutoff is needed only in the first emitted 
gluon because subsequent collinear emissions are already regulated by the angular ordering con- 
straint. 

Both approaches were compared when the rates for emission of a fixed number of resolved 
gluons, with a transverse momentum larger than a given resolution scale [ir, together with any 
number of unresolved ones, was performed in Ref. [8 |. The scale is constrained by the collinear 
cutoff and the hard scale, /i <C : <SC Q- It was found that jet rates in both multi-Regge (BFKL) 
and coherent (CCFM) schemes are the same. When coherence is introduced the singularities at 
uj = are stronger than in the BFKL approach but the extra logarithms cancel in the sum of all 
the contributions to the jet rates. Using a generating function for the jet multiplicity it is possible to 
prove that this is true to all orders in the coupling [J9][T0] . The same generating function is obtained 
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for BFKL and CCFM, with the p-th central moment of the jet multiplicity distribution being a 
polynomial in In -2. G f degree 2p — 1, showing that the distribution becomes relatively narrow in 
the limit of very small x and large Q/fiR \ 9 1 ■ 

In Refs. ifTTI the subject was developed further and subleading logarithms were included to cal- 
culate the minijet multiplicity associated to Higgs production at hadron colliders. In IPT21 it was also 
shown that for any sufficiently inclusive observable the CCFM formalism should lead to the same 
results as the BFKL equation. For the interested reader good reviews devoted to the implementa- 
tion of CCFM in Monte Carlo event generators can be found in, e.g., lfT4HT7l . For Monte Carlo 
implementations of the BFKL approach see Refs. (T3ljT6l lT8 - 22l. More recent results related to the 
implementation of models with unitarization in the CCFM formalism can be found in Refs. M231I24I 



3 Monte Carlo implementation of the CCFM evolution and nu- 
merical results 

The numerical implementation of the CCFM equation we use in the present work is the Monte Carlo 
event generator SMALLX developed by Marchesini and Webber in Ref. Il25l . In this code the CCFM 
gluon Green function, interpreted as the unintegrated gluon structure function J CCFM (k 0) k&,a;) 
when acting on the initial condition shown below, can be written in the iterative form 



f^ M (k a ,k b ,x) 



E 



6{x- x Q ) (k Q - k b ) A s (rj, q ) 6(1]- q ) 
0{q- z n l n ) A s (77, z n l n ) 5(x- x n ) <5 (2) (k a - k n ) . 



(8) 



As we mentioned in the previous section, q is the collinear cutoff for the first real emission. k 4 
corresponds to the transverse momentum of the exchanged gluons, with Xi being their longitudinal 
momentum fraction. l{ = qi/(l — zi) is the rescaled transverse momentum of the emitted gluons. r\ 
is the upper limit of the phase space integration. The Sudakov form factor reads 



As (U, Zi-ik-x) 



exp 



-2a, 



Zi — ili — 



dl 

T 



' dz 



1 - z 



where A is a cutoff for soft singularities. The final result is independent of A in the limit A 
probability of emission of each of the real gluons can be written in the form 



dPSi = A s (h, Zi-ik-i) P (zi,qi,ki) 9 (h - z^ik^^e 1 - 2* - 



A\ d 2 h 



L J nl; 



-dzi 



with the gluon splitting function being 

P(zi,qi,ki) 



A(z 2 ,qi,ki) 



1 



(9) 
0. The 

(10) 
(11) 



Our target is to compare this numerical implementation of the solution to the CCFM equation 
with the corresponding one in the BFKL approach. For this we have studied several distinctive 
distributions which will allow for the comparison. In this letter we focus on the leading order 
approximation in both approaches keeping the coupling fixed. We will investigate higher order 
corrections in future publications. 
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If we represent the CCFM solution in the iterative form of Eq. ([8]) we can test its convergence 
investigating what is the contribution to J CCFM stemming from a fixed number of real gluon emis- 
sions. For a given x = e~ Y when Y is larger configurations with more emissions have a bigger 
weight. This is seen in Fig. [T] where we also notice that the maximum of the distribution, with 
a s = 0.2, for Y = 2, 4, 6, 8 lies, respectively, at n = 4, 5, 7, 9. The chosen external momentum 




Figure 1 : Distribution in the contributions to the CCFM gluon Green function with a fixed number 
of emitted gluons, plotted for different values of the center-of-mass energy. 

scales are k a = 15, fcb = 20 GeV. In Fig. [2] we observe that the shape of the distribution is better 
reproduced by a Gaussian than by a Poissonian fit. 

/jv(8, 15 GeV, 20 GeV) 




5 10 15 20 



Figure 2: Gaussian (straight) and Poissonian (dashed) fits to the distribution in the contributions to 
the CCFM gluon Green function with a fixed number of emitted gluons, plotted for rapidity Y = 8 
and fixed values of the external momenta. 

As the gluon emissions take place it is more likely to populate regions of phase space far away 
from the transverse scales present in the initial condition. We have studied this point in Fig. [3] 
where we show the mean deviation from the central value of the typical transverse momentum of 
the emitted gluons from the external scales k a ,kb- We continue choosing a s = 0.2 and we now 
take k a = 15 GeV and kb = 20 GeV. The central lines show the maximum of the distribution in 
the momentum of the emitted gluon k for a given rapidity Y, showing all plots with different total 
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rapidities in a single figure normalizing them to the interval from zero to one. The two outermost 
lines correspond to Y = 8 with the total Y decreasing as we move closer to the central line. It is 



0.4 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 



y 



Figure 3: Distribution of transverse scales in the evolution with x of the CCFM equation. 

very interesting to see how the diffusion into the infrared region, where k 2 < k a kb, dominates the 
evolution. To quantify this statement we have calculated the ratio 



"^uv/ir — — (12) 



Area uv 
AreaiR 

of the area above (ultraviolet diffusion) over the area below (infrared diffusion) the central line. For 
Y = 2, 4, 6, 8 we obtain TZxjy/m = 0.42, 0.43, 0.44, 0.47. From these numbers we can see that the 
contribution of the infrared region is enhanced as we decrease the available energy in the scattering 
process. If we introduce a cut in the number of gluons considering only those contributions to the 
gluon Green function with more than 10 emissions then we obtain lZ\jy /jr = 0.39, 0.48, 0.51 for, 
respectively, Y = 4,6, 8. This indicates that high-multiplicity contributions are more dominated by 
infrared effects for low total rapidities than the low-multiplicity configurations, while converging to 
a more UV/IR symmetric structure for larger Y. This is natural since they dominate the full Green 
function in this region. 

When comparing with the BFKL results, it is interesting to investigate the collinear/anticollinear 
behaviour of the solution to the CCFM equation. This can be done by studying, with a fixed rapidity 
Y, the regions with large/small ratio of the external scales kb/k a . We have done this in Fig.|4]for 
different values of the 'reference scale" kb = 5, 10, 30 GeV. The main conclusion of this analysis 
is that at lower values of the external scale kb the gluon Green function becomes much flatter as a 
function of the variation of the other external scale k a . One is naturally tempted to relate this be- 
haviour to a possible approximate "conformal invariance" present in this limit. If we now introduce 
a cut in the number of emissions and keep only those contributions with more than ten gluons in the 
final state then we obtain Fig. [5] We can see that the approximate scale invariance present in the low 
kb cases remains in the low-multiplicity configurations, being a quite universal feature of the CCFM 
radiation. 

We now focus on the azimuthal angle dependence of our numerical solution. The azimuthal 
angle, 6, we are interested in is that formed by the two two-dimensional external vectors k a and kj,. 
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Figure 4: Collinear behaviour of the CCFM gluon Green function for different values of the refer- 
ence scale, kb. 
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Figure 5: Collinear behaviour of the CCFM gluon Green function for different values of the refer- 
ence scale, kb, considering the contributions of configurations with more than 10 emissions. 
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For this we extract the corresponding Fourier components using 

f% CFM (k a ,k b ,x) = jT ' ^-f CCFM (k a ,k b , x) co S (n9). (13) 

In Fig. [6] we observe that the n = component is the dominant one, with the n > components 
also growing with energy but at a slower pace. We will see that this is completely different to the 
BFKL case. If we impose the high-multiplicity cut then we can see in Fig. [7] that at larger Y the 
Green function is completely dominated by the many-gluon configurations. We can also conclude 
that the growth of all Fourier components is a common feature at any multiplicity. 
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Figure 6: Variation with rapidity of the different components of the Fourier expansion on the az- 
imuthal angle of the CCFM gluon Green function. 




Figure 7: Variation with rapidity of the different components of the Fourier expansion on the az- 
imuthal angle of the CCFM gluon Green function, considering the contributions of configurations 
with more than 10 emissions. 
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To wrap up this section we consider the full 8 dependence in Fig. [8] It is clear that the bulk of 
the configurations live in the region with k a "back-to-back" with kt, (6 = 0) with a Gaussian-like 
spread towards other angular settings. At high multiplicities the k a and k^ momenta tend to be more 




Figure 8: Full dependence on the azimuthal angle of the CCFM gluon Green function, 
decorrelated in the azimuthal angle as seen in Fig. [9] The combined (9, Y) dependence is shown 




Figure 9: Full dependence on the azimuthal angle of the CCFM gluon Green function, considering 
the contributions of configurations with more than 10 emissions. 

in Fig. [10] When only the high-multiplicity configurations are kept we obtain a slower variation in 
angles in Fig.fTTI 
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Figure 1 1 : The combined dependence of the CCFM gluon Green funtion on the azimuthal angle and 
rapidity, for fixed values of the modulus of the boundary momenta, considering the contributions of 
configurations with more than 10 emissions. 
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4 Monte Carlo implementation of the BFKL evolution and nu- 
merical results 

Our numerical solution of the BFKL equation is based on the following iterative representation for 
the LO gluon Green function: 



/ BFKL (k a ,k b ,Y) = e'M-^ ^)( ko - kb ) 



oo n 



•En a WnS*(*<- A2 ) / *e 2wS V'k-fe+^a, (i4) 

n=li=l J n i JO \ l=1 Jj 



where 



ui, 



o 




wo -(ki + ^k, -wo - ki+^k, . (15) 



\ \ i=i 

Here we have used this notation for the gluon Regge trajectory: 

^ = "i7j p (q -kF " -y l0g ^' (16) 

where t — —q 2 and A is an infrared regulator. The final result is independent of A in the limit of 
small A. The normalization of our numerical implementation corresponds to the following analytic 
form of the Green function: 

f BFKL (k k 6 Y)-^- T f^e» Y [^(^Y" e — (17) 

where 8 is the azimuthal angle between the k a and kb transverse momenta. In the case of elastic 
scattering n can be interpreted as a conformal spin in the unitary principal series representation of 
SL(2, C) (261. The eigenvalue of the BFKL kernel is 

X (n, 7 ) - 2*(l)-*( 7 +|)-*(l- 7 + ^). (18) 

The Monte Carlo analysis is useful because it allows for a more detailed study of different 
distributions. In particular, we can, as we did in the CCFM case, investigate the convergence of the 
iterative solution in terms of the contributions with a fixed number of emissions for a given value 
of Y ~ log 1 jx. We have shown this distribution in Fig. [T2]for Y = 2, 4, 6, 8, k a = 15 GeV and 
fcb = 20 GeV and a s = 0.2. We observe a very similar pattern to the one present in the CCFM 
case: a broadening towards larger number of emissions for large rapidities with the maxima being 
at n — 2, 5, 8, 11 for, respectively, Y = 2, 4, 6, 8. The main difference here is that in the BFKL case 
the distribution is clearly Poissonian. 

The diffusion picture in BFKL at leading order is well understood. We can average over the 
azimuthal angle: 

2tt 



/ BFKL (fc a ,fc b ,Y) = -L [ / BFKL (k a ,k b ,Y) 
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Figure 12: Distribution in the contributions to the BFKL gluon Green function with a fixed number 
of emitted gluons, plotted for different values of the center-of-mass energy. 



where x(l) has poles at integer values of 7 with the physical region corresponding to < 7 < 1. 
At asymptotic values of Y the dominant region of integration corresponds to 7 ~ i . The 7^1 — 7 
symmetry of \ indicates that there exists a symmetric diffusion into infrared and ultraviolet modes: 
since the pole at 7 ~ 0(1) corresponds to virtualities in the internal propagators smaller (larger) 
than the external scales, it leads to infrared (ultraviolet) diffusion. To be more precise: for very large 
values of a s Y we can write 

/™ (*.,*>,¥) = -^r/SfiV" 5 ^^ (20) 



irk a kb J 2iri \k 2 

and evaluate the integral using the saddle point approximation around the minimum of x 



\(-! - 41,.-2 + l-U,( 7 -i) +■■■ (21) 



to obtain the expression 



p FKL (k a ,k b ,Y) ~ _^ P a s 4io g 2Y 1 ( 22) 

with t = log (ka/k?). It is easy to verify that the function 

$(fc Q ,fc b ,Y) = k a k b p FKL (k a ,h,Y) (23) 
asymptotically fulfills the diffusion equation 

= 4 log2$ + 14(3^, (24) 



d(a s Y) ta sa dt 2 ' 
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which is independent of the sign of t and therefore, going back to the definition in Eq. (fTZt . 
T^-uv /ir = 1 in the BFKL case. It is interesting to note that this diffusive behaviour is only driven 
by the anomalous dimension 7 while the n ^ do not play any role. The fact that in the CCFM case 
the ratio Tvluv /IR tends to 1 for very large rapidities is in agreement with the fact in the small x limit 
the CCFM and BFKL approaches have the same asymptotic limit for the anomalous dimensions 
governing the scale variation of DIS structure functions. 

The collinear/anticollinear limits, at leading order, are not as interesting as in the CCFM case 
since the Green function has the functional form g{k a /k ) / (k a k ). Changing the reference scale fcj, 
does not bring any new features as we can see in Fig. Qj] where we chose Y = 6. When imposing a 
cut to only keep high-multiplicity contributions we obtain a similar result, see Fig. [T4l 

LO BFKL Green function - collinear behaviour 

°- 1 [ ' ' ' ' ' k b '=1QGeV j 

b 30 GeV : 





k a / k b 

Figure 13: Collinear behaviour of the BFKL gluon Green function for different values of the refer- 
ence scale k a 

Apart from the different diffusive behaviour, it is in the azimuthal angle dependence where we 
find more differences between the BFKL and the CCFM approaches. The expansion on Fourier 
components in the azimuthal angle can be written as 



/ BFKL (k Q ,k b ,Y) = £ /„ BFKL (fc a ,fc b ,Y)e^, 
with the coefficients being 



(25) 



n— — oc 



fn {ka, kb, Y) — 



chf 



■Kk a k b J 2ni \^fc2 



(26) 



When performing our Monte Carlo analysis we obtain these coefficients using the momentum space 
numerical solution projecting on angles, i.e. 



f^(k a ,k b ,Y) = £^/ BFKL (k a ,k b ,Y)cos(n ( 



(27) 
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LO BFKL Green function - collinear behaviour for large multiplicity 
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Figure 14: Collinear behaviour of the BFKL gluon Green function for different values of the refer- 
ence scale k a , considering the contributions of configurations with more than 10 emissions. 



In Fig. [15] we show how the convergence in n for the Green function is very fast. In this example 
we can see that 10 terms in the series are enough to reach a very good approximation to the final 
solution. This agrees with the general behaviour obtained with our Monte Carlo analysis in Fig. [16] 
The combined plots in rapidity and azimuthal angle for full and high multiplicities are given in 
Fig. [I7]and Fig. HU 

Very importantly, in the BFKL approach the Fourier projections in 8 have a very different be- 
haviour to those in the CCFM formalism. As we can observe in Fig.[l9]the n = projection also 
rises as in the CCFM case in Fig. [6] However, the non-zero n BFKL components do not rise with Y 
while the CCFM ones do rise, see, again, Fig. [6] This indicates that any observable sensitive to these 
higher angular components will have a completely different behaviour in both theories. There has 
been recent progress in the study of the azimuthal angle dependence within the BFKL approach, see, 
e.g., Refs. [27-31 36-38]. In future publications we will present observables where this different 
azimuthal angle dependence might help discriminate between BFKL and CCFM contributions. A 
natural candidate is the azimuthal angle decorrelation between jets produced in the central region of 
rapidity and a jet emitted along the direction of one of the hadrons in the Large Hadron Collider. 



5 Conclusions and scope 

In this letter we have compared two Monte Carlo implementations of the CCFM and BFKL 
formalisms for the description of small x observables. The main difference between them from 
the theoretical point of view is the introduction of QCD coherence effects in the CCFM equation. 
We have found that the symmetric diffusion into infrared and ultraviolet regions of phase space 
characteristic of the BFKL parton evolution is broken in the CCFM case, where the infrared scales 
play a dominant role. As our main result we have found that the higher Fourier components in 
the gluon Green function have a very different behaviour in both theories, rising with energy 
in the CCFM case and decreasing in the BFKL one. It will be very interesting to trace these 
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Figure 16: Full dependence on the azimuthal angle of the BFKL gluon Green function. 
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f(Y, 15 GeV, 20 GeV, <|>) 




Figure 17: The combined dependence of the BFKL gluon Green function on the azimuthal angle 
and rapidity, for fixed values of the modulus of the boundary momenta. 



f(Y, 15 GeV, 20 GeV, <|>) 




Figure 18: The combined dependence of the BFKL gluon Green function on the azimuthal angle and 
rapidity, for fixed values of the modulus of the boundary momenta, considering the contributions of 
configurations with more than 10 emissions. 
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LO BFKL Green function - Fourier components in azimuthal angle 
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Figure 19: Variation with rapidity of the different components of the Fourier expansion on the 
azimuthal angle of the BFKL gluon Green function. 

differences at an observable level [32-35] and to implement higher order corrections [36. 39—1471 to 
evaluate their effects on them. These lines of research will be the subject of our future investigations. 
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